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1 Introduction 

Infinite divisibility and their specialized versions, namely, self-decomposability and stability, 
have generated considerable interest among specialist in probability and statistics. There is 
huge literature devoted to studies of these topics. Books such as Loeve (1963), Linnik (1964), 
Feller (1966) and Lukacs (1970) have been instrumental in providing the audience with the 
basic material on these. More recent monographs such as Bondesson (1992), Sato (1999) and 
Steutel & van Harn (2004) have unified and studied further contributions to the expanding 
literature in this connection. 

In view of Kendall & Stuart (1963, Chapter 3), in which the relations between moments 
and cumulants are addressed in detail, it follows, under appropriate conditions, that the 
cumulants corresponding to infinitely divisible distributions that exist have some appealing 
features and have links with certain moments of Levy and Kolomogorov measures relative 
to these distributions. Ramachandran (1969, Theorem 8) has shown, in the univariate case, 
that for an infinitely divisible distribution the existence of the absolute moment of order 
a G (0, oo) is equivalent to the existence of its analogue for a certain truncated version of 
the corresponding Levy measure. (By a truncated version of a measure v on W (p > 1) 
we mean the restriction of v to some proper subset of W.) This result plays a crucial role 
in studies related to cumulants of infinitely divisible distributions, see, e.g., Gupta et al. 
(1994) and Steutel & van Harn (2004, Chapter IV, §7). Sato (1999, Theorem 25.3) has given 
a multivariate extension of Theorem 8 of Ramachandran (1969), involving the so-called g- 
moments. However, it appears that Theorem 25.3 of Sato (1999) is not tailored to meet the 
immediate needs for cumulant studies. 

In a recent expository article, Gupta et al. (2009) have unified the literature on infinitely 
divisible distributions with special reference to moments and cumulants. In the process of 
doing this, they have made several illuminating observations on the behavior of cumulants 
of univariate and multivariate infinitely divisible distributions, and have presented some 
new results in the area. Gupta et al. (2009) also poses an open problem on a multivariate 
extension of Theorem 8 of Ramachandran (1969). One of the main tasks of the present article 
is to deal with this problem; the problem that we have referred to here is of particular interest, 
especially if one is concerned with aspects of cumulants of multivariate infinitely divisible 
distributions. Interestingly, as a by-product of our solution to the problem, it follows that 
for any infinitely divisible distribution on W (p > 1), under a mild assumption, (in standard 
notation) the cumulant fc ri ,..., r exists if the moment fJ> ri ,...,r p exists (as a real number); in the 
univariate case, obviously, this result holds without requiring the distribution to be infinitely 
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divisible. 

In a somewhat different direction, there are questions relative to structural aspects of the 
multivariate hyperbolic distributions of Barndorff-Nielsen (1977) and their extensions with 
densities given by equation (7.3) in the cited reference; each of the distributions referred to 
here is indeed (in the notation of Barndorff-Nielsen (1977)) a mixture of N n (fi + u(3A,uA) 
with respect to u, where u follows a certain generalized inverse Gaussian distribution and 
/x, P and A are fixed with A nonsingular. (The extended versions have been termed the 
generalized hyperbolic distributions, especially in the univariate case by, e.g., Halgreen (1979, 
p. 14) and J0rgensen (1982, p. 37).) As claimed by Shanbhag & Sreehari (1979, p. 24), there 
exist members in the class of multivariate generalized hyperbolic distributions that are not 
self-decomposable. Specific examples illustrating that this is so can be found in, e.g., Pestana 
(1978, p. 54) and Rao & Shanbhag (2004, Example 3.3, Remark 3.4); Rao & Shanbhag (2004) 
consists of further information on the problem telling us, amongst other things, that there 
exist members also in the smaller class of multivariate hyperbolic distributions that are not 
self-decomposable. In this article, we attempt a comprehensive solution to a characterization 
problem that is linked with the question on the structural aspects of multivariate hyperbolic 
and multivariate generalized hyperbolic distributions being addressed. 

The paper is organized as follows. In Section 2, a generalization of Theorem 8 of 
Ramachandran (1969) to the case of multivariate distributions is provided in conjunction 
with several relevant observations and pertinent examples. In Section 3, a characterization 
theorem, based on the property of self-decomposability, is established for a certain class of 
mixtures of multivariate distributions, and its implications are emphasized. As mentioned 
before, the results of Section 2 are of importance in cumulant studies and the results of 
Section 3 throw further light on the structural aspects of multivariate generalized hyperbolic 
distributions and some related distributions. 

2 Criteria based on Levy measure for the existence of moments 
for multivariate infinitely divisible distributions 

From Levy (1954), or any other appropriate source such as Feller (1966, XVII. 11) or Sato 
(1999, Theorem 8.1), it follows that <fi is the characteristic function (ch.f.) of an infinitely 
divisible (i.d.) distribution on W 1 (p > 1) if and only if it is of the form 

m = exp ji < a,t > - l -Q{t) + £ ^<^> - 1 - ) » teW > ( 2J ) 
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where a is a real vector, Q is a nonnegative definite quadratic form, and v is a measure 
(referred to as Levy measure) on the Borel a-field of MP such that z/({0}) = and 

12 



\x\ 



1 + 11*1 



dv(x) < oo. 



It is easily seen that (2.2) is equivalent to the condition that 



(min{||x||, r}) 2 dvix) < oo for any fixed r G (0, oo) 



(2.2) 



(2.3) 



(Here, < •, ■ > and || • || denote respectively the usual inner product and usual norm on M p .) 
As observed by Gupta et al. (2009), using essentially the approach of Loeve (1963, 



Complement 9, p. 332), with r G (0, oo), we can rewrite (2.1) as 



(t) = exp ii < b,t > --Q(t) + / (e iKt ' x> -l-i<t,x>)du 2 (x)\ 



x exp <^ -z/i(M p ) + 



J<t,£C> 



di>i(x) 



t G 



(2.4) 



where z^i(-) = v ({x : \\x\\ > r} fl ■), i>2 = v — Ui, and (in obvious notation) 

|2 



a r + 



\x\ 



dv 2 (x) 



1 + \\x\ 



dui(x), r = l,2, ...,p. 



(Note that (2.2), or the equivalent condition (2.3), trivially implies that 6 r 's are well defined 
as well as that v is a a-finite measure and v\ is a finite measure.) 

We begin now by giving the following theorem that extends Theorem 8 of Ramachandran 
(1969) to the case of i.d. distributions on MP '. The theorem is clearly in a format which 
makes it easily applicable in cumulant studies, especially when the constants a r and (3 r , 
r = 1, 2, . . . ,p, defined in it are integers; the expressions for cumulants corresponding to an 
i.d. distribution on M p are related, under appropriate assumptions, to those for the moments 



of measure z/i referred to in (2.4), as observed, e.g., in Gupta et al. (2009). 



Theorem 1. Let X = (X\,X2,. . . ,X P ), p > 1, be a p-component i.d. random vector 



vr=l 



with ch.f. <p satisfying (2.1) (and hence also (2.4)), and let (3 r G [0, oo), r = 1, 2, . . . ,p. Then, 

E (^JJlA^I^ < oo for all a r G[0,/? r ], r = l,2,...,p, (2.5) 

if and only if 

Y\ \x r \ ar j dui{x) < oo for all a r G [0,/3 r ], r = l,2,...,p, (2.6) 
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vr=l 



where x = (x l5 x 2 , ■ ■ ■ , x p ) and V\ is as in (2.4). 



Proof. That (2.5) implies (2.6) is an obvious corollary to the first assertion of Theorem 



5.1 of Gupta et al. (2009). Now, to prove that (2.6) implies (2.5), we may proceed as follows. 



Since X is an i.d. random vector with ch.f. <p, in view of (2.4), we can see that <p is of the 
form 

<f){t) = <p 2 {t)exp{-X + X<f) 1 {t)}, t G W, (2.7) 



where A = u^W), 



hit) 



A_1 fnp ei<t ' X>d M x ), ^K", if A>0, 
1, teW, if A = 



and 



b 2 (t) = exp <j i < b,t > + 



/ ( e i<t, x > _ i _ i < t, x >) dv 2 {x) \ , te 



Essentially, adopting the approach of Feller (1966, p. 534) or others such as Sato (1999, 
Lemmas 25.6, 25.7) or Steutel & van Harn (2004, Chapter IV, §7) given in the univariate 
case, it can be easily seen that the distribution corresponding to the ch.f. <j) 2 has a (full) 



moment sequence and, hence, satisfies the analogue of (2.5). Since this disposes of the case 



of A = trivially, assume then that A > 0, where A is as in (2.7). From (2.7), we see, in this 
case, that 



(t) 



3=0 



t G 



(2.8) 



which, in turn, implies (on noting, in particular, that is a mixture of ch.f.'s (p 2 (t)((f)i(t))^ , 
t G W, for j G {0, 1, . . .}, where the mixing distribution is Poisson with mean A) in view of 
Fubini's theorem that, for each a r G [0, /3 r ], r = 1, 2, . . . ,p, 



E ( jj \X r 



vr=l 



3=0 



v r=l 



2r 



(fe) 
It 



k=l 



(2.9) 



where (X 21 , X 22 , . . . , X 2p ) and (xg } , xg\ . . . , X™), k 



1,2, ... ,j, are mutually 



independent random vectors with the first one having ch.f. <ft 2 and each of the remaining 



ones having ch.f. <pi. (In (2.9) we adopt the convention that the summations with respect 
to k equal if j = 0.) 



Assume now that (2.6) holds. In view of (2.9) and (2.6), it easily follows that, for each 
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3=0 
oo 



(fc) 
Ir 



3=0 

oo 



,r=l k=0 
j P 



,fc=0 r=l 



Y (k) 



where 



max 



3=0 



e n(i 

\r=l 



\r=l 



(2.10) 




< OO, 



and, for notational convenience, we denote X 2r by X$ for each r. (We appeal here to, 
amongst other things, the relevant independence of the vectors concerned and the elementary 
inequality, for a G [0, oo), y k G R, k = 0, 1, . . . , j, that | Y,l=o Vk\ a < (J + 1 )" max < fc <i |2/fc| Q < 



(j + 1)" XlLo bfcl"' an d also view the last summation with respect to j in (2.10) as the 



expected value of a function of a Poisson random variable.) Obviously, (2.10) yields the 
validity of (2.5). This shows that (2.6) implies (2.5). Hence, Theorem 1 follows. □ 



The following corollary of Theorem 1 is essentially a version of Theorem 8 of 
Ramachandran (1969); the cited author considers v\ with a \x\ > 1" in place of u \x\ > r". 

Corollary 1. Let X be an i.d. random variable with ch.f. 0, and let (3 G [0, oo). Then, 



if and only if 



E(|X|0) < oo 



/ \x\P dui(x) < oo, 
Jul 



where and v\ are as in (2.4) corresponding to p = 1. 



Proof. The corollary follows readily from Theorem 1 by taking p = 1 because if 
H is a finite measure on the Borel a-field of K then j R \x\@ dp,(x) < oo if and only if 
L \x\ a dp,(x) < oo for all a E [0,(3]. □ 

Remark 1. We point out that Theorem 1 is not covered by Theorem 25.3 of Sato (1999) 
on g-moments (relative to the case of r = 1); that this statement is true is obvious from the 
example (i.e., Example 4) given in Remark 5.5 of Gupta et al. (2009). However, to illustrate 
that neither the function J 



r=l 



\x r \ ar , met in the statement of Theorem 1, nor its modified 
version given by the function Y\fr=i \ x A ar ^-{A), x G MP, r > 0, with A = {x : \\x\\ > r} and 
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1(A) as its indicator function, is assured to be submultiplicative (where the terminology refers 
to that of Definition 25.2 of Sato (1999, p. 159)), we may consider the following example: 

Example 1. Let p = 2, x = (x\, x 2 ) G ^ 2 and r > 0, and consider 
gi(x) = \xiX2\I(A) and #2 (as) = l^i^l, x G M 2 . 
For x, y G IR 2 such that Xi > r, 2/2 > r, £2 = .x^ 1 and yi = y^ , we have 
5f r (a? +y) = \(xi + yi)(x 2 + y 2 )\ > x x y 2 , r = 1, 2, 

and 

r (a;) = p r (y) = 1, r = l,2. 
Consequently, it is impossible to have here for each r G {1,2} a constant a T > such that 

g r (x + y) < a r g r (x)g r (y), r = l,2, 

for all x,y £ IR 2 . This supports the claim that we have made above. (It is also now clear 
that the functions max{|xiX2|, c} with c > are not submultiplicative.) 

Remark 2. We may modify the example given in Remark 5.5 of Gupta et al. (2009) 
to shed further light on aspects of Theorem 1. In particular, we can show, with appropriate 
modifications to the example referred to, that if one or more of certain /3 r 's are positive then 



(2.6) with "[0, p r ] (in respective places)" replaced by "(0,/? r ]" does not imply (2.5) with the 
same change, and demonstrate some curious phenomena of (Xi, X 2 ) in this connection. That 
this is so, is evident from the information supplied by the following two examples: 

Example 2. Let (Xi,X 2 ) be an i.d. random vector with ch.f. such that 
</>(t) = exp{-A + A^(t)}, t = (t u t 2 ) G M 2 , 
where A G (0, 00) and if) is the ch.f. of a random vector (Yi, Y 2 ) satisfying 

(Y 1 ,Y 2 ) ^ (V^- 1 ), 

where V is the modulus of a standard Cauchy random variable (and "=" denotes the equality 
in distribution). Clearly, we have then for a\,a 2 G [0, 1] such that < |«i — a 2 \ < 1, i.e., 
whenever cui, a 2 both lie in (0, 1] (or both lie in [0, 1)), 

/ \x 1 \ ai \x 2 \ aa du 1 (x)< J x^xfdvix) = \E(V lai ~ a21 ) < 00. 

Jm 2 J(0,oo) 2 
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However, in this case, we have 



E(\X 1 X 2 \)=E(X 1 X 2 ) = oo and E (\X r \) — E (X r ) — oo, r = l,2, 

violating the condition that E |X 2 | Q2 ) < oo for all ai,a 2 G (0,1]. (Appealing to 

Theorem 1, we can also see in this case that E (|ATi| Ql |X 2 |° 2 ) < oo for each an, a 2 G [0, 1).) 

Example 3. Let {X\,X-i) be as in Example 2 with the exception that ip in this case 
refers to the ch.f. of a random vector (Yi, Y 2 ) satisfying 

(Y U Y 2 ) ^ (V^V- 6 ), 

with 7 G (0, 1), 5 G (0, oo) and as a positive stable random variable with left extremity 
(i.e., inf{x : F v (x) > 0}, where F Vl denotes the distribution function (d.f.) of V 7 ) zero and 
characteristic exponent 7. It now follows that, for a±,a 2 G [0, 1] with a 2 7^ 0, 

/ \xi\ ai \x 2 \ a2 dv^x) < [ x^xf dv{x) = AE (V^~ a2&) ) < 00. 

(That the expectation appearing above is finite is seen, e.g., from Bondesson (1992, p. 85) 
or from Steutel & van Harn (2004, p. 246).) Also, it is clear that, in this case, X\ and X 2 
are nonnegative random variables such that E(Xi) = ~E(X\X 2 ) = 00, with X 2 possessing a 
(full) moment sequence. 

The following theorem enables us to understand the mechanism of Theorem 1 better; 
the theorem addresses, amongst other things, the problem posed in Remark 5.6 of Gupta et 
al. (2009), and its proof that we have produced here is adapted partially from the proof of 
Theorem 5.1 of Gupta et al. (2009). 

Theorem 2. Let X = {X\, X 2 , . . . , X p ) and /3 r , r = 1, 2, . . . , p, be as in Theorem 1, but 
for that there is an additional restriction now that F(X r = 0) < 1, r = 1, 2, . . . ,p. Then, 



(|2.5|) is equivalent to the condition that 

v 



e ni^i^ < °° • ( 2 - n ) 



K r=l 



Moreover, we now have 



F S t[\ x r\ > \ > 0. (2.12) 



Proof. Trivially, the theorem is true for p — 1. We shall follow the method of induction 
with respect to p to prove it, noting that each subvector of X is i.d. Assume then that p > 1 
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and the theorem holds in the case when p' , with p' G {1,2, ...,p — 1}, appears in place of 



p in (2.11) and (2.12). The theorem is clearly valid if X is p-variate normal and hence it 



is sufficient if we show that the theorem holds when v is non-null. Consequently, we can 
assume, without loss of generality, that ^(R^) > 0, where E + = [0,oo), and observe that 



the theorem follows if, under the assumption, it is proved just that (2.12) is valid and, in 



view of Theorem 1 and the prevailing symmetry, that (2.11) implies (2.6) with W replaced 
by R p + . 



Let A = z/l(R+) > 0; then (2.4) implies that, for n — 1,2, .. ., we have p-variate ch.f.'s 



a 2 and 03 such that 



t) = ln Mt)^'(t) + (i - 7n) t e » p 



(2.13) 



where 



(*)) 



n-l 



exp{-A + A0i(t)} 



t G M p , 
, t G M p , 



A 



-i 



A<t,x> 



duAx), t g W, 



with 7„ = e~ x X n /n\. (Note that A and 0i considered here are not implied to be the same as 
those met in the proof of Theorem 1.) 

Letting (X n , Xu, • • • , Xi p ) and (X^i \ X^\ • • • , X^), n = 1, 2, . . ., denote independent 
random vectors with ch.f.'s 0i and 02™^? n = 1,2, . . ., respectively, we get from (2.13) that 

(2.14) 



E ( H \X r f" ) > ln E ( H \x lr + Xg> 



\r=l 



vr=l 



n = 1,2, 



and (noting especially that, for each n, 02 x ^(t) = 1 (t)0 2 ri ' ) (t), t G 



p nix r i>o > 7n p mix£ +i; 



,r=l 



,r=l 



> 



n = 1,2, 



(2.15) 



If k < p of the Xi r , r = 1,2, ...,p, are equal to zero almost surely and the remainder 
satisfy the condition that P{X lr > 0} > 0, then we can take without loss of generality that 
P{X n = X 12 = . . . = X lk = 0} = 1 and P{X lr > 0} > 0, r = Jfe + l, k + 2, . . . ,p; we take here 
k = to mean that P{Xi r > 0} > for all r G {1, 2, . . . ,p}. For a sufficiently large integer 



no, the distribution corresponding to 2 nas a ^ least one support point, say (c\,C2, ■ ■ ■ ,c p ) 
with c r 7^ for each r G {1, 2, . . . ,p} and c r > for each r G {& + 1, k + 2, . . . ,p}; to see this 
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note that, by assumption, 



V{ Ill^l > \ > if k > 1, 



(2.16) 



and F{X n = X n = . . . = X lk = 0} = 1 together with P{X ir > 0} > 0, r = jfe+1, k+2, ...,p. 
(Obviously, (2.16) is a consequence of the assumption in the inductive argument, especially 



because we have now that, for each n, (X&\x$, . . . , X$) = (X u X 2 , . . . , X k ).) 

In view of the aforementioned observation on the existence of the support point 
(ci, C2, • • • , c p ) of the distribution relative to 02™°^ with the stated properties, it follows that 



P{lj;» ) 6A r ,r = l,2 ) .., J ,}>0, 



(2.17) 



where 



J 1 . 



(fsf) if c r <0, 
(f,^) if c r >0. 



(2.18) 



Consequently, by (2.15) we get that (2.12) is valid, and, by (2.14) in conjunction with (2.11) 



that, for some constant rj £ (0, oo) and A^s as in (2.18) 



oo > E J] \X r fr > V E [H \X lr + Xi: 0) " I Xt 0) £ A r , r = 1, 2, . . . ,p . (2.19) 



vr=l 



vr=l 



c;|'" | < OO, 



v.r=l 



where c* = |c r |/2, r = 1, 2, . . . ,p. Since, for a r £ [0, /3 r ], r = 1, 2, . . . ,p, 

\X lr \ ar < (4)°^ \X lr + 4f r , r = 1, 2, . . . ,p, 



In view of the properties of (X n ,X 12 , . . . , X lp ) and (X^ 0j , X^ oj , . . . , Xg oj ), p. 19) implies 
then that 



(2.20) 



it is hence obvious from (2.20) that 



E ( Y[ \X Xr \ ar I < oo for all a r £ [0, (3 r ], r = 1, 2, . . . ,p, 



vr=l 



asserting that, as required for the completion of the proof of the theorem, (2.6) with in 
place of MP, is met. Hence, Theorem 2 follows. □ 

Remark 3 (i). Theorem 2 is not valid if the assumption of i.d. is dropped. This is 
obvious from the following examples in which we use an indirect approach based on the first 
assertion of Theorem 2 to ascertain that the X considered are non-i.d.: 
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Example 4A. Let V and V* be independent random variables such that E(|V|) = oo 
and V* is {0, l}-valued Bernoulli. Define X 1 = VV* and X 2 = V(l - V*). Clearly, the 
random vector X = (X 1 ,X 2 ) is such that P{X r = 0} < 1, r = 1,2, with P{XiX 2 = 0} = 1 
and hence with E(|XiX 2 |) = < oo. Also, in this case, obviously E(|Xl|) = E(|X 2 |) = oo. 
That X is non-i.d. and the claim of Remark 3 is valid follows then trivially from the first 
assertion of Theorem 2. However, it is interesting to observe that in this example, we have 
E(|Xi| ai |X 2 | aa ) = < oo for all a u a 2 G (0,1]. 

Example 4B. Let 7 G (0,1) and X = (Xi,X 2 ) be a random vector such that X\ is 
a positive random variable with E(X 7 ) = 00 and X 2 = Xf 1 almost surely. Then, clearly, 
we have E(|X 1 ||X 2 |) = K(XiX 2 ) = 1 < 00. However, in this case, it is not even true 
that EdXiHXal" 2 ) < 00 for all a u a 2 G (0,1), since E(|X 1 | 7+<5 |X 2 |' 5 ) = E(Xj) = 00 if 
5 G (0, 1 — 7). That X is non-i.d. and the claim of Remark 3 is valid follows again trivially 
from the first assertion of Theorem 2. 

Remark 3 (ii). That the X vectors dealt with in Examples 4A and 4B are non-i.d. 
can also be shown via alternative approaches without involving the findings of Theorem 2; 
in the remainder of this remark, and in Remark 3 (iii), we illustrate as to why this is so. Let 
(U, W) be a 2-component random vector with nonnegative and nondegenerate components 
U and W such that P{£/H^ = c} = 1 for some (nonnegative) constant c. Since, in this 
simple corollary to Theorem 2 of Shanbhag (1988), it follows that the distribution 
of (U, (UW) 1 / 2 , W) is indecomposable, it is obvious then that the distribution of (U,W) is 
indecomposable; to see this, note, in particular, that (UW) 1 ^ 2 is degenerate. Consequently, 
we have the distribution of X in Example 4A in the case when V is nonnegative and that 
of X in Example 4B to be indecomposable and hence non-i.d. 

Remark 3 (iii). Let (U, W) be as in Remark 3 (ii), but for a modification that U and 
W in this case are not necessarily nonnegative and also that c is allowed here to be negative. 
Applying essentially a simpler version of the argument used in Shanbhag (1988) to prove its 
Theorem 2, one can see that the distribution of (U, W) is decomposable if and only if for 
some 61, b 2 ^ with b\ — Acb^ 1 > 0, and some a, f3 G (0, 1), 



a(3 



if 



x 



(a 1 ,a 2 b 1 b 2 ), 



P{(U,W) =x}= < 



(1- 



a) /3 



if 



if 



x 



x 



(a 2 ,a 1 b 1 6 2 ), 
(-02, -aifef 1 ^), 



I (1 



a)(l 



0) if 



x 



(-ai, -a 2 &i 1 b 2 ), 



where 
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(reducing, when c = 0, to a\ = b± and a 2 = 0, respectively). (Note that the "if" part 
of the assertion follows easily since, under the relevant conditions, there exist independent 
2-component random vectors Y {1) and Y {2) so that F{Y {1) = (a 1} a 2 b^b 2 )} = a, F{Y (1) = 
(a 2 , ai b^b 2 )} = l-a, F{Y® = (0, 0)} = p, F{Y^ = -b 2 )} = 1-/3 and Y^+Y® is 
distributed as (U,W).) Clearly, the characterization met here implies that the distribution 
of (U, W) is indecomposable if U and W are nonnegative, a result referred to in Remark 
3 (ii), and also that the distributions of X appearing in both Examples 4A and 4B are 
indeed indecomposable and hence non-i.d; to see the validity of the claim concerning the 
distributions of X, note, in particular, that each of these examples has E(|Xi|) = oo with 
X meeting the assumptions relative to (U,W). 



, X p ), p > 1, with d.f. 



Remark 4. For any p-component random vector X = (X\,X 2 , 
F, irrespectively of whether or not it is i.d., (2.12) is met if and only if the support of F 

c„) so that Ylr=i \ c r\ > 0- Suppose now that, for some positive 



includes a point (ci,c 2 , 
integer j, 



5 ^p) 



X 



5> 

n=0 



(n) 



where Y {n \ n = 0,1,... , j are independent p-component random vectors such that, for each 
n, given any support point (ci,c 2 , . . . , c p ) of the d.f. of Y^ n \ there exists a support point 
(di,d 2 , . . . ,d p ) of the d.f. of Y^n'(^n)=o ^ n ^> f° r which Yll=i \dr\ and c r d r > for each 
r with c r 7^ 0. In this case, Ylr=i \ c r + d r \ > 0, from which we can see that (2.12) is met, 
and, following essentially the relevant steps in the proof of Theorem 2, appearing (2.17) 
onwards, we can further see that (2.11) implies (2.5) with 7" in place of X respectively 
for n = 0,1, ... ,j, and hence, in view of the inequality referred to under brackets below 



(2.10), in the proof of Theorem 1, that the first assertion of Theorem 2 holds; this is so 



irrespectively of whether or not X is i.d. (To see, especially, that if (2.11) holds for this 
X, then, for each n, Y^ n ' satisfies (2.5), it is sufficient, by symmetry, taking, in obvious 
notation, without loss of generality that P{l^ (n) G ({0}) fc x (0, oo) p ~ k } > 0, to check that 
this is so for a vector with distribution the same as the conditional distribution of Y^ given 
that Y {n) e ({0}) fc x (0,oo)^ fc .) Consequently, we are led now to a new version of Theorem 
2. 

The following corollary, which is obvious from Theorem 1 and 2, presents an extended 
version of the first assertion (i.e., the crucial assertion) of Theorem 5.1 of Gupta et al. (2009). 

Corollary 2. Let X = (X\,X 2 , . . . , X p ) and (3 r , r = 1, 2, . . . ,p, be as in Theorems 1 and 
2. Then, (2.11) is equivalent to (2.6). 
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Remark 5. In statistical literature, whenever the covariance between two random 
variables is mentioned, the relevant random variables are often assumed to be square- 
integrable. The following simple example, which indeed is a slight variation of Example 
2 met above, illustrates that the covariance may be well-defined even when this standard 
assumption does not hold and the joint distribution of the random variables is i.d.: 

Example 5. Let (Xi,X 2 ) be an i.d. random vector with ch.f. <p such that 

0(t) = exp{-A + \ip{t)}, t = t 2 ) £ R 2 , 

where A £ (0, oo) and ip is the ch.f. of a random vector (V, V" 1 ) with V as a positive random 
variable such that its square equals the modulus of a standard Cauchy random variable. 
Then, we have 

E(|XiX 2 |) = E(XiX 2 ) < oo and E (|X r |) = E (X r ) < oo, r = l,2. 

However, in this case, it is obvious that E (Xf) = E (Xf ) = oo, hence we have the validity of 
our claim. (A closer scrutiny of the Example 5 above tells us, actually, in view of Theorems 
1 and 2, something more. Indeed, the random variables X% and X 2 that we have considered 
in this example satisfy the condition that E (|Xi| ai |X 2 |° 2 ) = E(X 1 ai X 2 ° 2 ) < oo for each 
cci, ct2 £ (0,2].) Incidentally, it may be worth noting here that, in view of the corollary to 
Theorem 2 of Shanbhag (1988), implied in Remark 3 (ii) above, any random vector relative 
to ch.f. if> of the present example gives us an example of a random vector for which the 
conclusion of Remark 5 is valid with "i.d." replaced by "indecomposable". 



3 A characteristic property relative to multivariate self- 
decomposability 

Following Urbanik (1969, p. 92), let us define a distribution on W (p > 1) to be self- 
decomposable (s.d.) if the corresponding ch.f. <p satisfies the condition that, for each 
c£ (0,1), 

<j>{t) = (j){ct)(j) c {t), teW, (3.1) 

where <ft c is a ch.f. Essentially, from Sections 8 and 11 of Chapter XVII of Feller (1966), it is 



then clear that if (3.1) is met then both <p and <p c are i.d. ch.f.'s; the partial information that 
(f) and 4> c are nonvanishing also follows, in effect, from the argument given in Lukacs (1970, 
p. 162) to show that this is so in the univariate case. (As a by-product of this, it follows 
that if is nonvanishing or, in particular, i.d., then is a s.d. ch.f. if and only if </>(£) j '(f>(ct) , 
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t G MP, is an i.d. ch.f. for each c G (0, 1).) If a distribution on M. p is s.d., we refer to the 
corresponding ch.f., or a random vector with this distribution, also as s.d. 

As pointed out in Section 1, the concept of generalized hyperbolic distributions is 
well documented, especially in the univariate case. Shanbhag & Sreehari (1979, p. 24) 
have claimed that there exist members in the class of multivariate generalized hyperbolic 
distributions with /3 ^ (where the distribution referred to is absolutely continuous with 
probability density function (p.d.f.) given by (7.3) of Barndorff-Nielsen (1977)) that are not 
s.d. The original example implied in the cited paper to illustrate this remained unpublished 
though its version was later revisited in Pestana (1978, p. 54). More recently, Rao & 
Shanbhag (2004, Example 3.3, Remark 3.4) have produced a simpler example in support 
of the claim and have made some further relevant comments on the issue. There is also 
Example 6.8 in Gupta et al. (2006) that enlightens one with certain related information. 

Before discussing the results of the present section, we may give two further pieces of 
information that are of relevance to these: 

(i) Although Feller (1966, XII. 8) is concerned with properties of s.d. distributions in the 
univariate case, that the properties of <j) and </> c referred to above are met in the multivariate 
case is essentially implied by Feller (1966, XII. 11); also, it is of interest to note here that, 



if cf> satisfies (3.1), then every linear combination of the components of the corresponding 



random vector has its ch.f. satisfying the univariate version of (3.1), implying immediately 
that and C (of the multivariate case) are nonvanishing. 

(ii) To be more precise, in the notation of Section 1, we refer in this paper to mixtures 
of N n (n + u/3A,uA) as generalized hyperbolic if the mixing distribution relative to u is 
generalized inverse Gaussian, i.e., if it has p.d.f. of the form 

f(u | A,£,V0 = C(A,e,V)w A - 1 exp{-(^- 1 + <M/2}, u > 0, 

with £, ip > for which max{£, ip} > 0, C as the normalizing constant depending on the 
modified Bessel function of the third kind, and A so that / °° w A_1 exp{ — (£w _1 + ipu) /2}du < 
oo. Obviously, the class of hyperbolic distributions is a subclass of the class of these 
distributions, see, e.g., Barndorff-Nielsen (1977). 

The following theorem and its corollary subsume some of the major observations that are 
made by the examples cited above. 

Theorem 3. Let V be a positive random variable such that 

E (e sV ) = exp ( / (e sv - 1) v~ {a+l) g(v)dv\ , Re(s) < 0, (3.2) 

U(0,oo) J 
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with a G [0, 1) and g as a bounded decreasing nonnegative real function on (0, oo) satisfying 



v g[v) 



dv < oo. 



'(0,oo) l+« 

Also, let X = (Xi,X 2 , ■ ■ ■ , X p _i), p > 2, be a (p — l)-component random vector independent 
of V, with X r , r — 1, 2, . . . ,p — 1, as independent (nondegenerate) symmetric stable random 
variables with characteristic exponents 7 r G [1, 2], r = 1, 2, . . . ,p— 1, respectively. Then, the 
p-component random vector 



Z= V,VnX 1 ,...,Vp-iX. 



p-1 



is s.d. if and only if 



a 



p-i 

r=l 



(3.3) 



(3.4) 



Proof. In view of (3.2), it follows that the ch.f. of Z given by (3.3) is of the form 



(Pit) = exp |y (j^-vZr=l VltH-iP- _ , t = (t 1? t 2 , . . . ,t p ) G M p , 

(0,OO) (3.5) 
with A r > for all r G {1, 2, . . . ,p — 1}; this is clear since the ch.f. of each symmetric stable 
random variable X r , with characteristic exponent j r , is of the form <p r (t) = exp {— A r |t| 7r }, 
t G R, and we have 



p-i 

(t) = E ( e i,ir ]jE 



t G 



From (3.2), it is obvious that <fi is i.d. Also, on appealing to Fubini's theorem, it is now clear 



that the Levy measure in the present case is concentrated on (0, 00) x ]R P_1 and is absolutely 
continuous with respect to Lebesgue measure on M? with Radon-Nikodym derivative h such 
that the restriction to (0, 00) x ]R P_1 of h is given by 



p-i 



h(y) = yl {a+1) g(yi) U 



r=l 



y=(yi,V2,...,y P ) G (0,oo)xM p - 1 , (3.6) 



where f r denotes the p.d.f. of X r and is implied by the inversion theorem to be bounded 
(continuous). Recalling then that <fi is s.d. if and only if, for each c G (0,1), </>(£) / 4>(ct) , 
t G MP, is i.d., we can claim that Z is s.d. if and only if 



h{y) > c p h(y/c), y G (0, 00) x 



np-l 



(3.7) 
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for all c G (0,1). Since each symmetric stable distribution is unimodal with vertex (see, 



e.g., Lemma 5.10.1 of Lukacs (1970)), (3.6) then implies that (3.7) is met for the required c 
and y if and only if (3.4) is valid; observe in particular that if (3.4) is not valid, then (3.7) 
is violated for y with y r+ i = o(yl^ r ), r = 1,2, ... ,p — 1, and yi sufficiently small. Hence, 
Theorem 3 follows. □ 



Remark 6. In view of (3.2) directly, or, trivially, as a corollary to Theorem 3, it 
follows that the random variable V considered in Theorem 3 is s.d. Positive stable (with 
left extremity zero), gamma and inverse Gaussian random variables provide us with some 
specialized versions of V met here; in these cases, we have g(v) oc e~ with A = for stable 
and A > otherwise, and, also, have the parameter a respectively as positive, equal to 0, 
and equal to 1/2. Also, in view of the closure property (under weak convergence) of the class 
of s.d. distributions on W, as a corollary to Theorem 3, it now follows that the "if" part of 
the theorem referred to remains valid without the assumption that g be bounded. 

Corollary 3. Let V and X r , r = 1, 2, ... ,p — 1, be as in Theorem 3, and, additionally, 
suppose that we have p > 3 and a — p + 2 + Y7r=tlr l < 0. Then, there exist real c r , 
r — 1, 2, . . . ,p — 1, such that the (p — l)-component random vector W given by 

W = (c x V + V%X 1} c^V + V^X P ^ 

is not s.d. 

Proof. Clearly, it follows via a standard argument that the sequence of the random 
vectors {W n : n — 1, 2, . . .}, where, for each n > 1, 

w n = (-V + V*X U . . . , -V + v^x p _ 2 , V + -V^Xp^) , 
\n n n J 

converges in probability and hence in distribution to the random vector W* given by 

W* = (v^Xx, V^X p _ 2 , V 

Since, by Theorem 3, we have that W* is not s.d., appealing to the closure property 
(under weak convergence) of the class of s.d. distributions on M. p ~ l (p > 3), we can 
readily conclude that, for some n > 1, W n is not s.d. This, in turn, implies that the 
assertion of the corollary is true; note that if, for some n, W n is not s.d., then so also is 
(ly + VnX 1} . . . , \V + V^X p _ 2 , nV + V^X p ^). Hence, Corollary 3 follows. □ 

Remark 7. Although, the problem of finding non-s.d. multivariate hyperbolic and 
multivariate generalized hyperbolic distributions of Barndorff-Nielsen (1977), touched upon 
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in this article, prompted us to establish Theorem 3, that this latter theorem is of interest 
in its own right, is obvious. However, it may be worth emphasizing here that Corollary 3, 
which is a corollary to Theorem 3, identifies certain members of the class of multivariate 
generalized hyperbolic distributions of Barndorff-Nielsen (1977), that are non-s.d.; this 
follows on noting especially that the specialized version of the corollary in the case of 
a G {0,1/2}, g(v) = exp{— Xv}, v > 0, and 71 = 72 = ■■■J p -i = 2, concerns these 
distributions. The examples of non-s.d. distributions given by Corollary 3 are obviously in 
the spirit of those discussed earlier in Pestana (1978, p. 54) and Rao & Shanbhag (2004, 
Example 3.3 & Remark 3.4, pp. 2882-2883). 

Corollary 4. Given an integer p > 2, there exists a p-component random vector of the 



form of (3.3) that is not s.d. such that all its lower dimensional subvectors are s.d. 

Proof. Given p > 2, choose, e.g., 7 G (1, 2] and a G [0, 1) such that a = (p — 2)(1 — 7 _1 ) 
and, hence, satisfying also that a — (p — 1)(1 — 7 _1 ) < 0. Then, Theorem 3 implies that in 
the special case of 7 r = 7, r = 1, 2, . . . ,p — 1, with a as stated, the random vector Z is not 
s.d., but, for each r G {2, 3, . . . ,p — 1}, the (p — l)-component subvector of Z that does not 
include the r-th component of Z is s.d. Also, since the specialized version in this case of Z 



with its first component deleted satisfies (3.1) trivially for each c G (0, 1), it is obvious that 
this is s.d. Hence we have the corollary. □ 

Corollary 5. If a = 0, the random vector Z is s.d. if and only if 71 = 72 = ... = 7 p _i = 
1, i.e., if and only if Xi, X 2 , . . . , X p _i are Cauchy (up to scale changes) random variables. 

Proof. The result follows trivially from Theorem 3; also, the "if" part of the assertion 



is immediate on noting that the concerned ch.f. satisfies (3.1) for each c G (0, 1). □ 

Remark 8. Clearly, any random vector Z* is i.d. if its ch.f. is of the form 

p-i 



I iti — A r |t r+ i| 7r j , t — (t±,t2, . . . ,t p ) 
V r=1 / 



with, p > 1, a G [0, 1), X r > and 7 r G (0, 2] (i.e., p, a, X r and 7 r are as in Theorem 3 but 
for that 7 r is now allowed to lie in (0, 1)), and t>* as the function defined by the left hand side 
of (3.2); note that Z* — Z, where Z is as in (3.3) if 7 r G [1,2], r = 1,2, . . . ,p— 1. Moreover, 
essentially as in Example 6.8 of Gupta et al. (2006), it is now seen that each projection of 



Z* is s.d. if 71 = 72 = ... = 7 P -i = 2 and, additionally, g in (3.2) is completely monotone, 
or, equivalently, by Bernstein's theorem, denotes the Laplace transform of a measure on the 
Borel cr-field of IR + ; this follows since in the present case also we have the distribution of 
V to be a member of the generalized gamma convolution family of Bondesson (1992), on 
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observing that if g is a Laplace transform as above, then so also is v~ a g(v), v G (0, oo). 
However, by Theorem 3, we have, in this case, obviously Z* (and hence Z) to be non-s.d., 
unless p = 2 and a G [1/2, 1). 

Remark 9. The example provided in the proof of Corollary 4 gives us the existence 
of yet another i.d. distribution on W p (p > 2) that is not s.d., of which each projection 
is s.d. This, in conjunction with the information provided by Corollary 4 and Remark 8, 
compares well with certain findings of Levy (1948) and Shanbhag (1975); the results in the 
cited references show us, amongst other things, that there exist indecomposable distributions 
on IR P (p > 2) with all marginals (univariate or otherwise) and projections as i.d. Further 
material of relevance to the findings referred to here has appeared or was cited in, e.g., 
Davidson (1973), Kendall (1973), Shanbhag (1974, 1976, 1988), Ostrovskii (1986) and Letac 
(1992). 

Remark 10. Since the class of distributions on MP that are s.d. is closed under 
convolution, from an observation in Remark 3.4 of Rao & Shanbhag (2004, p. 2883), it 
is clear that any ch.f. of the form 

< (iciti + ic 2 t 2 - Ait? - X 2 tl) , t=(h, t 2 ) G M 2 , 

with t>* as in Remark 8, (01,02) ^ (0,0) and Ai, A2 > 0, is s.d., provided that a G [1/2,1) 
and g is completely monotone; to see this, use the fact that, in this case, v ^^ 1 ^ 2 ^g(v ), 
v G (0, 00), is completely monotone. The cited remark of Rao & Shanbhag (2004) also gives 
some further relevant information on the subject. 



Remark 11. In view of Sato (1999, Example 25.10, pp. 162-164), by ( |OD , it follows by 
Fubini's theorem (in the notation of Theorem 3) that if a — p + 1 + Ylr=i TV" 1 < 0' then, for 
each (5 G (0,p— 1 — a — YllZilr 1 ) an d G (a /(a + (3), 1), we have a function G on (0, 00) 
such that for all x G (0, 00), 



G{x) = [ yj^^^f J] \y r+1 \Adv{y) < 00, (3.8) 
J (.art \r=i / 



where v is the Levy measure relative to the distribution of Z (with Z as in (3.3)). For each 
c G (0, 1), in obvious notation, the analogue of G with respect to ch.f. 4>(ct), t G MP, 



can easily be seen to be given by c («+/ 3 -!'+ 1 +Er = i 7 r ^ e G(x/c), x G (0, 00), with G as in (3.8); 
obviously, since, for each c G (0, 1) and x > 0, G^(x) > G(x), it then follows that, for none 
of c G (0, 1), we have, in this case, <p(t)/<p(ct), t G MP, to be i.d. Consequently, we have now 
an alternative argument for proving the "only if" part of Theorem 3. 
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Remark 12. Suppose, we define (in the notation in (3.2)) 



6* = sup \e > : / v e - a - l g{v)dv < oo 1 



'[l,oo) 

Then, by Corollary 1, and the result from Sato (1999, Example 25.10, pp. 162-164), we have, 
in view of Fubini's theorem, that, for each (a%, a 2 , ■ ■ ■ , a p ), p > 2, such that a r G (—1, 7*_i), 
r = 2, 3, . . . ,p, and «i + £]£ =2 a r 77-i € [0, 0*), 



/ p 

Eizf 1 n\z 

^ r=2 



< OO, 



where Z = (Z\, Z 2 , . . . , Z p ), p > 2, as in (3.3) but for a modification that, in this case, j r , 
r — 1, 2, . . . ,p — 1, p > 2, are allowed to be less than 1, and 7*, r = 1, 2, ... ,p — 1, are so 
that 7* = 7 r if 7 r < 2 and 7* = 00 if 7 r = 2. In this remark, as in Remark 11, we have come 
across arguments essentially in the spirit of those met in Section 2 of the paper. 

Remark 13. Theorem 3 does not hold if the assumption that 7 r G [1,2], r = 
1, 2, . . . , p — 1, is dropped. This is obvious from the following example: 

Example 6. Let {<ft n : n — 1, 2, . . .} be a sequence of i.d. ch.f's on IR 3 , such that 

Mt) = < (iti-t 2 2 -^\t 3 py t = (ti,t 2 ,t 3 ) et 3 , 

with t>* as in Remark 8, a G [0,1/2) and 7 G (0,2/(3 - 2a )). Clearly, the sequence 
{4> n } converges to a ch.f. on IR 3 , which, by Theorem 3, is not s.d.; note that the limiting 
random vector (Y~i,Y 2 ,Y 3 ) in this case is so that the ch.f. of (Y 1 ,Y 2 ) equals v*(iti — t|), 
(ti,t 2 ) G M 2 , which is indeed non-s.d. by Theorem 3. Consequently, appealing to the closure 
property (under weak convergence) of the class of s.d. distributions on IR 3 , we can claim 
that there exist, for some large n, ch.f.'s (f) n that are not s.d. on IR 3 despite the fact that 
a -2 + 2- 1 + 7 - 1 > 0. 

Remark 14. Applying a result of Zolotarev (1954) (which has appeared also as Theorem 
5.8.4 of Lukacs (1970)), in conjunction with a standard result stating that each stable 
distribution is unimodal (an obvious consequence of a result of Yamazato (1978) stating 
that each s.d. distribution is unimodal), one can easily see, with appropriate scrutiny of 
Levy measures, that there exist i.d. ch.f.'s on IR 2 of the form 



<f>{t)=v*Jit 1 -\ 1 \t 2 y), t=(t 1 ,t 2 )E 



d2 



with v* as in Remark 8, Ai > and 7 G [1/2, 1), that are not s.d. for certain a and g 
(e.g., if a < 1 — 7 and g is identically equal to a constant). This sheds further light on the 
observation of Remark 13. 
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